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Hysteretic resistance spikes in quantum Hall ferromagnets without domains 
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We use spin-density-functional theory to study recently reported hysteretic magnetoresistance p^x 
spikes in Mn-based 2D electron gases [Jaroszyriski et al. Phys. Rev. Lett. 89, 266802 (2002)]. We 
find hysteresis loops in our calculated Landau fan diagrams and total energies signaling quantum- 
Hall-ferromagnet phase transitions. Spin-dependent exchange-correlation effects are crucial to sta- 
bilize the relevant magnetic phases arising from distinct symmetry-broken excited- and ground-state 
solutions of the Kohn-Sham equations. Besides hysteretic spikes in p^x, we predict hysteretic dips 
in the Hall resistance pxy Our theory, without domain walls, satisfactorily explains the recent data. 
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Two-dimensional electron gases (2DEGs) under strong 
magnetic fields exhibit fascinating physical phenomena; 
the mostly widely known of these being the integer and 
fractional quantum Hall effects [IJ. Spontaneous mag- 
netic order in quantum-Hall systems is yet another inter- 
esting possibility. Quantum Hall ferromagnetism arises 
from the interplay of the Zeeman, Coulomb, and thermal 
energies within the macroscopically degenerate Landau 
levels of the 2DEG 0, Q ■ Landau level crossings Q offer 
a convenient means to probe symmetry-broken quantum- 
Hall ferromagnetic transitions. At crossings, opposite- 
spin levels can benefit from Coulomb exchange to form 
spin-ordered states at low temperatures 0] ■ 

Many groups have investigated quantum Hall ferro- 
magnetism in the integer and fractional quantum-Hall 
regimes by inducing Landau level crossings via tilted 
magnetic fields, density and level tuning via gate elec- 
trodes, hydrostatic pressure 0], and the s-d exchange- 
induced level bowing in Mn-based 2DEGs d, [3| . These 
studies d, @, 0] find ubiquitous "anomalous" peaks and 
hysteretic spikes in Shubnikov-de-Haas measurements of 
the magnetoresistance pxx- The hysteretic behavior here 
means that the spikes appear at distinct magnetic fields 
as the field is swept up and down. In the fractional 
quantum-Hall regime recent experiments show that these 
features follow from the hyperfine coupling between elec- 
trons and nuclei In the integer quantum Hall regime, 
the hysteretic spikes have been suggested to arise from 
charge transport along long domain-wall loops acting as 
one-dimensional ("percolating") channels in the 2DEG 
near an Ising-like quantum-Hall ferromagnet transition 
Q. Though appealing, this description is largely quali- 
tative: no magnetotransport quantities (e.g., pxx) have 
been calculated so far accounting for domain walls. 

Alienate picture for the hysteretic spikes. In this Let- 
ter we apply the Spin Density Functional Theory (SDFT) 
[lot implemented via the Kohn-Sham (KS) scheme in a 
Local Spin Density Approximation (LSDA) to determine 
the electronic structure of a 2DEG, which we then use 
in a linear response model to explicitly calculate pxx 
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FIG. 1: Hysteretic magnetoresistance pxx, Hall resistance pxy, 
and total energy E, as a function of B for a quantum well 
[inset in (a)]. For the up^down B sweeps, hysteretic spikes 
and dips appear in pxx and pxy, respectively, at B ~ 5.8 T (a). 
For increasing temperatures the hysteretic loops become less 
pronounced {pxy not shown) and disappear above a critical 
temperature Tc ~ 2.1 K (b). The total energy E (c) shows 
hysteretic loops [see blow-up of the loops in (d)] because the 
KS equations have two distinct solutions (ground and excited 
states), e.g., for < B < _B". These are quantum-Hall 
ferromagnet phases with distinct spin polarizations (Fig. 2). 
The curves in (b) are displaced vertically for clarity. 



and Pxy- For concreteness, we focus on the ex per iment 
of Jaroszyhski et al. in Mn-based 2DEGs In- 
terestingly, we find hysteretic behavior in pxx and pxy, 
Figs. l(a)-(b), without taking into account domain walls. 
This behavior follows from our system having two self- 
consistent KS solutions ~ for a same set of parameters 
- with distinct total energies as shown in Figs. 1(c)- (d), 
see e.g. range < B < B^. As we discuss later on, 
these solutions correspond to ground and excited states 
describing distinct quantum Hall ferromagnetic phases of 
our interacting 2DEG, and, most importantly, provide a 
concrete example for a theorem of Perdew and Levy on 
the existence of excited states from ground-state density 
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functionals [l3|. These phases comprise differing sets of 
conducting states contributing to the magnetotransport 
and hence have distinctive Pxx^^ and PxyS. The hysteresis 
then arises because the 2DEG can go through a different 
sequence of magnetic phases (i.e., the system can become 
trapped into distinct local minima) as the B field is swept 
up and then down, e.g., phases I to II in the up sweep 
at i? ~ B" and phases III to IV in the down sweep at 
B ~ B'i # B", Fi g. 1 (d). We predict hysteretic dips (and 
peaks [131) ill Pxy \1M Pig- 1(a)] and that the spikes shift 
to opposite directions (Fig. 3) in samples with positive 
and negative g factors, for increasing tilt angles of B. 

Mn-based system. We consider a CdTe quantum well 
between Cdo.8Mgo.2Te barriers, with three evenly spaced 
Cdi_2;pMna;pTe monolayers ("Mn barriers") in the well 
region Q, Fig. 1(a) (inset); Xp is the planar concentra- 
tion of Mn. Adjacent to the barriers lie two symmet- 
ric n-doped regions [l^. In an external field B, the s-d 
exchange interaction between the electrons in the well 
(2DEG) and those of the localized d orbitals of the Mn 
gives rise to a spin-dependent electron potential 
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where Noa is the s-d exchange constant, -B5/2 is the 
spin-5/2 Brillouin function, x{z) and Tq are the effec- 
tive Mn profile and temperature respectively, and 
±1 (or t, 4) denotes the electron spin components. 
The structural confining potential of the well is assumed 
square Vw{z) = vo['d{—L/2 — z)-\-'d{L/2-\-z)] with depth 
vq; Q{z) is the Heaviside function. The structural poten- 
tial of the Mn barriers is Vb{z) = vix{z); vi is the barrier 
height and x{z) the nominal gaussian Mn profile. 

Kohn-Sham approach. We use the SDFT [l3| formu- 
lated in the context of the effective-mass approximation 
of semiconductors. Within the finite-temperature formu- 
lation of Mermin 17, 1^, we obtain the KS equations 
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where m is the effective mass, i — 1,2... the band index, 
and v'^ff{z] [n|,nj) the effective single-particle potential 
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In Eq. ^ Vs{z) = Vw{z) -f Vh{z), Vh{z;[n]) is the 
Hartree potential, calculated by solving Poisson's equa- 
tion, and Vxc{z; [n|,nj) is th e sp in- dependent exchange- 
correlation (XC) potential [19|. The motion in the 
xy plane is quantized into Landau levels with ener- 
gies En = {n l/2)huJc, n = 0,1,2,... and u)c = 
eB/m [i^]. The total wavefunction is i^in kyi-'^^V^ ^) ~ 
^ exp[ikyy)ifin{x)Xi^ {z), where (finix) is the n-th 



harmonic oscillator eigenfunction centered at xq = 
— hky/mujc and ky is the electron wave number along 
the y axis; Ly is a normalizing length. This de- 
coupling of the z and xy motions follows from the 
uniformness of the total electron density within the 
xy plane: because each electron can be anywhere 
within the plane, we use the average total electron 
density n{z) = E.,n,ky,.JluL7L; 1 1 ^kkf^xdy = 

J2un,ky,a, \x1'iz)\'^fi,n in Poissou's equation, instead of 
n(x,y,z) [here f^^ is the Fermi fimction and is a 
normalizing length]. This procedure makes the 2DEG 
uniform thus rendering a separable KS set. 
We assume that the KS eigenvalues 
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where gp,BBaz/2 is the ordinary Zeeman term [g: effec- 
tive Lande factor), describe the actual electronic struc- 
ture of our 2DEG. This assumption is, in principle, un- 
justified within DFT: the individual KS orbitals repre- 
sent states of a fictitious non-interacting electron gas 
in an effective potential, Eqs. ^ and ([3]). With this 
assumption, however, we satisfactorily explain observed 
hysteretic phenomena in 2DEGs 0, [l2l ■ 

System parameters. In our simulations we use (see 
Ref. 0] ): m/irio = 0.099, dielctric constant e/eo = 
10, g = —1.67, A^oa = 220 meV, quantum-well width 
L = 100 A, spacer width Ls = 200 A, n2D = 2.97 x lO" 
cm^^, number of Mn monolayers Ni, ~ 3, Xp ~ 0.115, 
vo = 248.1 meV, vi = 1183. 5meV, To = 0.47 K, and 
assume a diffusion length £ ^ 4.67 A for the Mn profile. 

Figures 1 and 2 show our theoretical results for p^x , Pxy 
l(a)-(b); total energy E l(c)-(d). Landau level fan dia- 
gram 2(a), spin-resolved electron densities nl^'^ 2(b), and 
spin-polarization ^ 2(c). Remarkably, all of these quanti- 
ties show abrupt hysteretic changes near 5.8 T as the B 
field is swept up and down. The discontinuous nature of 
these features follows from the 2DEG undergoing (quan- 
tum) phase transitions in which its degree of spin polar- 
ization C suddenly changes, 2(b)-(c), e.g., in the down 
sweep the 2DEG becomes highly spin-polarized near 5.8 
T [2(c)]: a quantum Hall ferromagnetic phase transition 
takes place. Our calculated fan diagram 2(a) and total 
energy 1(d) corroborate this scenario: the opposite-spin 
levels e\ Q, e{ i suddenly cross near p, thus lowering the 
total energy E [see III to IV discontinuity in 1 (d)] . This 
is similar to the Giuliani-Quinn instability where the 1*^'- 
ordcr transition is due to the gain of exchange energy 
in the ferromagnetic state Q]- Our calculation also in- 
cludes correlation which somewhat reduces the exchange 
effects. The role of the s-d exchange [Eq. ((1])] in our 
system is to cause "level bowing" (B < 2 T) thus induc- 
ing opposite-spin Landau level crossings, see e.g. cross- 
ings for B < 4 T in 2(a). When crossings occur near 
p, spin-dependent XC effects may induce phase transi- 
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FIG. 2: Landau-level fan diagram [Eq. |4] (a), spin-dependent 
electron densities tIj^ (b), and spin polarization (c) of 
our 2DEG [Fig. 1(a)]. In (a) we plot the center of each 
Landau level with a phenomenological gaussian broadening 
r = 0.36\/BmeVT"^/^ [U. The s-d induced non-linear be- 
havior of the energy spectrum {B < 2T) allows for non-trivial 
Landau level crossings near the chemical potential fi (curve 
at ~ 74meV). Inset (a): blow-up of the hysteretic crossings 
between the states with n = 1 and spin down {e{ j) and that 

with n = and spin up {el q) at B ~ 5.8 T for the up^^down B 
sweeps. The discontinuous change in the spin-polarization (" 
signals a transition between distinct quantum-Hall ferromag- 
net phases. We have adjusted Xp and F so the levels would 
cross at ~ 5.8 T as in the experimental data 0. 



tions (e.g. at 5.8 T); we do not find phase transitions 
in a Hartree calculation. The phases here are those of 
an Ising-like itinerant 2D ferromagnct with Pauli mag- 
netization (n-f — ni)iiB = n2DCt^B- Next wc show how 
hysteresis arises in our SDFT calculation. 

Excited and ground Kohn-Sham states. Using the 
constrained-search definition of the ground state func- 
tional Ey [n] (defined on domain of densities constructed 
from any wave function), Perdew and Levy [f^l have 
shown that ". . . every extremum of Ey [n] represents the 
density ni{r) and the energy Ei of a stationary state. 
The absolute minima represent the ground states, and 
the extrema lying above the minimum represent a subset 
of the excited states. " These authors have also proved 
that some of the self-consistent solutions of the KS equa- 
tions extrcmize E^, [n] , provided that these solutions obey 
ground-state Fermi statistics [aufbau principle) with a 
single chemical potential - this is a necessary and suffi- 
cient condition. These theorems have been generalized 
to SDFT [2lj and hence hold in our system. 

Our simulations show that for a certain window of 
magnetic field [e.g., < 5 < in Fig. f(d)] the 
nations ^ have indeed two self-consistent solu- 
23 with distinct total energies and spin polariza- 
18| - both satisfying the necessary and sufficient 
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condition above. At each B in this range, one of these 
two solutions is stable (true minimum, ground state) and 
the other metastable (local minimum, excited state). In 
practice, these two states (phases) are separated by an 
energy barrier psj which may trap the system in the 
metastable states during the up^down B sweeps thus 
giving rise to hysteretic loops, Figs, f and 2. 

Magneto-transport. We obtain the longitudi- 
nal and Hall resistances from the conductivity 
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a ), calculated within the self- 



consistent Born-approximation model of Ando and 
Uemura [nj. For short-range scattercrs, a^x = 
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Ilj. We model the extended Landau states by a 
gaussian 5oxt(£) of width Foxt = 0.25 meV [24 1. 

Hysteretic resistance spikes and dips. The ordinary 
Pxx peaks in Fig. f (a), e.g., at 2.5, 3.1, and 5 T, are due 
to subsequent single Landau levels crossing the chemical 
potential p and appear between plateaus in p^y. Though 
involving the crossing of two opposite-spin levels at /i, the 
Pxx spike at B ~ 5.8 T has the same physical origin as 
the ordinary peaks: electrons in partially filled crossing 
levels are susceptible to scattering ( "dissipation" ) which 
increases pxx- This spike, however, corresponds to a dip 
in a plateau of pxy- as the two Landau levels cross near /x, 
the number of conducting channels Uc "fluctuates" thus 
making pxy dip (Aric > ) [or peak (Aric < 0)] toward an 
adjacent plateau consistent with n^. Whether pxy dips or 
peaks and whether the spikes and dips are stronger in the 
up or down sweeps depend on the details of the crossings 
at p, Fig. 2(a) - these features are, however, hysteretic as 
Uc differs in the up and down B sweeps. For increasing 
tilt angles 0, Fig. 3, the hysteretic spike shifts to lower 
fields while the ordinary Shubnikov-de-Haas maxima do 
not. Figure 3 uniquely identifies ordinary (single Landau 
levels crossing p) and "anomalous peaks" (two Landau 
levels crossing p) in pxx- 

Hysteresis & critical temperature. Figure 4(a) shows 




60/50/40/30/20/0 ^^^^2.97 x 1 0^ ^ cm'^ 

T = 0.33K 




6 7 

Bcos(9) (T) 



FIG. 3: Magnetoresistance pxx vs B cos{9) for several tilt an- 
gles 6 between the B field and the growth direction. Similarly 
to the data in 01, the spike shifts to lower fields as d increases, 
while the ordinary Shubnikov-de-Haas peaks do not. 
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FIG. 4: Temperature dependence of the pxx spike positions 
Be (a) and amplitudes (b) for up^^down B sweeps. Inset 
in (a): ABc = Be ~ versus T from which we extract a 
critical temperature Tc\ for T > Tc the hysteresis disappears. 
The asymmetric shape of our hysteretic loops [Figs. l(b)-(d)] 
is clearly manifest in (b). 



the hysteresis in the pxx spike position Be as a function 
of T: _B" > for the up ^down B sweeps. By plot- 
ting ABc = i?" — i?c versus T we can extract a critical 
temperature above which AB^ = 0. We find [inset in 
4(a)] Tc = 2.1 K which is comparable to the experimen- 
tal value 1.3 K [3]. The amplitude of the spike is also 
hysteretic, Fig. 4(b). Here, however, only the up-swcep 
behavior agrees with the data 0, HI] . 

Further comparison with experiments. We also repro- 
duce the non-hysteretic p^x peak at ~ 3.2 T (single level 
crossing n) seen in Q. We find that the peak at ~ 2.8 
T arises from two levels crossing /i and predict that 
it shifts upward (as opposed to Fig. 3 here) as in- 
creases because g > in while g < in 0]. However, 
we do not find any hysteretic behavior here. Essentially, 
the non-integer filling factors near the opposite-spin level 
crossings in [0] (as opposed to 0]) result in exchange en- 
ergy gains not enough to induce phase transitions. 

Final remarks. (i) Quenched-disorder-induced do- 
mains - not domains arising in metastable states [ p. 12511 
- can lead to transport anisotropics in 2DEGs [26|, l27| . 
(ii) The relevance of the hyperfine coupling to the hys- 
teretic phenomena has been recognized experimentally 
only in the fractional quantum Hall regime. We account 
for disorder effects only via the broadening of the Landau 
levels and neglect the hyperfine coupling and domains al- 
together. Our successful description of hysteretic (quan- 
tum Hall) ferromagnetic phenomena in 2DEGs highlights 
the power of DFT in a non-conventional application. 
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